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Abstract 

For the pedestrian observer, financial markets look completely random with erratic and 
uncontrollable behavior. To a large extend, this is correct. At first approximation the dif- 
ference between real price changes and the random walk model is too small to be detected 
using traditional time series analysis. However, we show in the following that this differ- 
ence between real financial time series and random walks, as small as it is, is detectable 
using modern statistical multivariate analysis, with several triggers encoded in trading sys- 
tems. This kind of analysis are based on methods widely used in nuclear physics, with large 
samples of data and advanced statistical inference. Considering the movements of the Euro 
future contract at high frequency, we show that a part of the non-random content of this 
series can be inferred, namely the trend-following content depending on volatility ranges. 
Of course, this is not a general proof of statistical inference, as we focus on one particular 
example and the generality of the process can not be claimed. Therefore, we produce other 
examples on a completely different markets, largely uncorrected to the Euro future, namely 
the DAX and Cacao future contracts. The same procedure is followed using a trading sys- 
tem, based on the same ingredients. We show that similar results can be obtained and we 
conclude that this is an evidence that some invariants, as encoded in our system, have been 
identified. They provide a kind of quantification of the non-random content of the financial 
markets explored over a 10 years period of time. 



1 Introduction 



The random walk model of price changes in financial time series has been so durable because 
it is nearly correct. At first approximation the difference between real price changes and the 
random walk model is too small to be detected using traditional time series analysis [Q]-[3]|. 
More precisely, when looking at large samples of data, some features appear that break the 
random walk approximation. For example, the statistics of price distribution at small time 
scales is not Gaussian but governed by non-extensive statistics BUO. We can also detect large 
range correlation in the absolute returns, which mean that persistent behaviors exist that are 
not embedded in the random walk model (TJ-Ol, which can be seen as a consequence of the 
non-extensive statistics J6]|. 

More explicitly, the non-extensive formalism provides an expression for the probability den- 
sity function P q of price returns x at a given time scale r [|2l[3l : 



P 9 (x, r) = {1 - P{t) [(1 -q)(x- x(r)) 2 ] } + (1) 



where q is a real parameter representing the degree of non-extensivity (q — > 1 in the Gaussian 
limit), Z q {r) is a normalization constant and f3{j) is a scale parameter. Also, is propor- 
tional to the variance of the distribution. In the expression of P q , the subindex + indicates that 
P q (x, r) = if the expression inside the brackets is non-positive. In general, for real markets, 
using large samples of data, the q index can be found to range from 1.2 to 1.7. This repre- 
sents intuitively the degree of the resulting anomalous diffusion from the underlying interaction 
among financial trades. Under certain approximations, regarding a free diffusion process, we 
can write (HH: 

0(t) -2 

Super diffusion occurs for q > 1. 

In Fig. [H we illustrate this behavior on the return distributions of the Euro future contract, 
sampled in 5 minutes units. The similarity in shape for two different years, 2002 where the 
Euro contract was moving up and 2005 where it was globally down is observed. These shapes 
correspond effectively to q > 1. However, any deviations from the Gaussian limit can not 
be detected at a local level, when we observe the market in a short window of time. Our first 
statement is still valid. At first order of observation, prices in financial markets behave randomly 
and it remains impossible to predict whether the next price movement will be up or down. 

In the following, we show that this difference between real markets and random walks, as 
small as it is, is detectable using modern statistical analysis with hypothesis testing, even when 
we observe the market locally. In particular, it is detectable once we wan build a trading systems 
on the basis of multivariate analysis and hypothesis testing I17UI41. Indeed, tools of statistical 
physics have been proven to be efficient in many areas, like extracting the average properties 
of a macroscopic system from its microscopic dynamics, even if approximately known. The 
same holds for financial systems. Even though it is difficult or almost impossible to write down 
the microscopic equation of motion that drives prices at each instant, it is possible to extract a 
relevant statistical information, that makes sense to take decisions at a local level. 
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In a first part, we exemplify this issue on the behavior with time of the Euro future contract 
(EC) at high frequency. We show that we can infer the non-random content of the EC erratic 
behavior using a multivariate analysis embedded in a trading algorithm. 

In a second and third parts, we examine different markets, largely uncorrelated to the Euro 
future, namely the DAX and Cacao future contracts (labeled as FDAX and CC). The same 
procedure is followed using a trading system, based on the same ingredients. For the system 
running on the CC, we use the same system as built on the Euro future. We show that similar 
and good results can be obtained on this variety of markets and we conclude that this is an 
evidence that some invariants, as encoded in our system, have been identified. 

2 Strategy on the Euro future 
2.1 Data sets and Data treatment 

We use five minutes sampling of the EC time series, from January 2000 till August 201 1, which 
makes 839k quotes that we use to build the trading system. We conserve only the close of 
each quote. This large sample of data points is necessary to infer statistical properties with a 
high confidence level, as shown in the following. Also, in the context of this analysis, the fine 
tuning of the time series with a five minutes resolution is useful to focus on possible intermittent 
behavior of the series at small scales (five minutes), that could disappear at larger scales. 

A typical quote of the EC is like 1.3802. The unit of the last digit is what we call a basis 
point. For example, we consider that a price movement from 1.3802 to 1.3805 corresponds to a 
price change of 3 basis points. More precisely, if we buy the contract at time Tl (on the quote 
Ql) at 1.3802 and sell this contract at time T2>T1 (on the quote Q2) at 1.3805, then this trade 
corresponds to a gain of 3 basis points (without fees). To keep the procedure as close to reality 
as possible we consider fees of two times the slippage, which means that this trade is counted 
in our approach as a trade of 1 basis point (net of fees). 

A fundamental issue in the analysis is to break the data samples in three parts, that we call 
in-sample, out-sample and live-sample. The decomposition is done as follows: 

(i) 2000-2007: in-sample 

(ii) 2008-2009: out-sample 

(iii) 2010-2011: live-sample 

What is the interest of this decomposition of the data series? The idea is that we intend to build 
a trading system on this series. This means that we intent to design an algorithm that will take 
decisions like buy or sell 1 EC contract at a given quote. This decision at a given quote will be 
based on multivariate analysis, as mentioned in the beginning. In order to process this way, we 
need a data sample on which the algorithm is built and all parameters of the algorithm are fitted. 
Therefore, this data sample needs to be large in order to be relevant statistically. This sample is 
called in-sample (i) and is defined as the period 2000-2007. 
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The second sample, called out-sample (ii), is used as a validation stage. All algorithms built 
on (i) are expected obviously to give satisfactory results on (i). However, as parameters of the 
model are fitted on the sample (i), there is no guarantee that the model could behave properly 
on another data sample. If it does so, this means that the algorithm is not a pure artifact and 
contains a part of the real dynamics of the market. This is the purpose of the sample (ii), defined 
as the period 2008-2009. If the trading system built on (i) fails on (ii), it is rejected and another 
algorithm is designed. Note that we have other intermediate validation stage to make the full 
process more robust: we come back on this point later in the article. Also, note that there is no 
guarantee at this level that what we describe in this paragraph is possible. 

Finally, once we have obtained an algorithm that works on (ii) and satisfies our robustness 
tests, if any, we observe it on what we call the live-sample (iii), defined as the period 2010- 
201 1 . Our building process is made to guarantee at this step the good functioning of the trading 
system and that's what we show in the following. 

Note that if we can drive the analysis to this last step and if it works, it is a clear proof of 
our claim of the previous part on a specific example (EC): the difference between real markets 
and random walks, as small as it is, is detectable using modern statistical analysis in multivari- 
ate analysis. The multivariate approach refers to the number of parameters introduced in the 
definition of triggers for trades decisions along the EC series. 

2.2 Strategy Reconstruction 

The basic elements of the algorithm design can be much simple. The gross featurse of a trend- 
following strategy are exposed in lfT5l . Let us note P(t') the value of the price series at time t' . 
An exponential moving average 4>(t) of length memory r can be defined on this time series, as: 



From being initially with no position, a trend-following system buys one share when </> 
reaches a given value $ and stays long until <\> hits the value — $, at which point the system 
sells back and takes the opposite position, and so on. A complexity can easily be added to this 
mechanism by defining intermediate thresholds to break positions taken by the system [fT5l . 
The trade distribution for this simple theoretical system is given in Fig. |2l Obviously, there is 
no possibility with such a simple algorithm to reconstruct a profitable strategy over ten years of 
high frequency data. 

However, we can use the ground idea of this mechanism, namely trend-following, in build- 
ing a more complex architecture. We use four different memory lengths and consider crossing 
of this exponential moving averages, as potential triggers for trade decisions. Not all moving 
averages are used for each decision. The choice is based on ranges of volatilities. Indeed, we 
have observed that there are some transition domains in volatilities of prices where it is prefer- 
able not to trade or to branch more stringent triggers. An important idea in our structure is also 
some exit conditions based on extreme conditions in profit, either on positive or negative. These 




(3) 
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last conditions depend also on the time window on which the non nominal cumulative profit is 
realized. 

On a fundamental ground, our system architecture is just a refinement of a basic trend- 
following strategy. In addition, the system has learned how to play with volatility ranges to 
trigger decisions and how to protect over-profits realized for example in high volatile periods. 
Therefore, if we can show that this strategy leads to profitable results (net of fees), it will be 
a proof of the validity of the trend-following hypothesis on the market, taking into account 
multivariate tests to activate the trend follower. 

Let us note that the use of moving averages is a powerful experimental method to access to 
the non-trivial statistical texture of a time series. If we consider 2 standard moving averages of 
lengths Ti and T 2 > T x , with AT = (T 2 - T\)/T\, then the density of crossing points of the 2 
averages is given by: 

p ^_L [Ar( l_ AT )]"-i (4 ) 

where H is the Hurst exponent, that characterizes the persistence or anti-persistent of the data 
series ifToll . 

We have 8 parameters optimized on the in-sample (i). The optimization is performed in 
order to achieve the best Sharpe ratio. Results are shown in Fig. [3] We present the behavior 
with time of the EC contract itself as well as the cumulative equity of the designed trading 
system (expressed in basis points, net of fees). We observe the nice behavior of the equity, 
increasing with time, which shows that the strategy is profitable and coherent with respect to 
different market regimes. The bottom plot in Fig. [3] corresponds also to the running of the 
trading system, but this time on the randomized in-sample. Exactly, we have added to each 
quote value of the data series (i) a random number that ranges between —10 and +10 times the 
slippage of the EC contract. And we run the trading system on this series, which leads to the 
bottom plot of Fig. [3] This randomization is necessary as we do not want the trading system 
to be dependent on the point-to-point correlation and also the model must be flexible to absorb 
distortion of the data series. This is what we observe in Fig. |3] (bottom): the system is robust 
against randomization of the data series. The degradation observed on the overall profit is not 
dramatic and the equity is still much reasonable. Note that all systems designed that have failed 
at this stage have been rejected. 

Before considering the out-sample stage, we have an intermediate essential step of valida- 
tion of the trading system. To ensure that the system is robust, we need more that the random- 
ization of the data series. We need to distort the strategy itself in many ways: for example, force 
an exit of given trade at a given time, do not execute randomly some trades, delay the execu- 
tion of orders by several quotes, execute an order but at a wrong price, with a prejudice for the 
trading system, multiply the fees (slippage) by a factor 2,3 or 4 etc. Thus, we have a list of 128 
stress tests and for each case, the trading system is run and a result is obtained. All this must be 
done on the original data series (i) and on its randomized version. In all cases, we must observe 
that the system is stable and robust. This is shown in Fig. HI where we present the Sharpe ratio 
for all 128 stress tests considered. We do not provide the equity in each configuration. We 
summarize each case by one entry in Fig. |4l as a value of the Sharpe ratio for the case under 
study. The idea is that the robustness is ensured if we do not observe pathological values in the 
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Sharpe ratios, even for the more extreme stress tests. This is what we observe in Fig. |U with an 
average value of 2.2 and a RMS of 0.8. In all configurations, the model stays reasonable. 

By this method, we have also shown that the trading system does not depend on the fine 
tuning of any of the fitted parameters. Otherwise, a few stress tests would have failed deeply. 
On the contrary, our strategy depends weakly on any of its inputs, which gives a lot of flexibility 
on all variables of the system with always a profitable result obtained. 

At this step, it is not unreasonable to claim that we have designed a robust algorithm. How- 
ever, a new validation stage is determinant using the out-sample (ii). This is a decisive test as we 
are running on new data, that the system does not know, in the sense that parameters have been 
fitted on another set of data. In principle parameters are robust as we have already explored 
many configurations for the data series and the system. However, the out-sample test will kill 
all systems that still have some elements of over-fit in their construction. Indeed, such systems 
fail to give good results when running on the sample (ii) and are rejected. This is what happens 
for most of the systems that can be designed if the input ideas are not carrying decisive features 
of the inside dynamics of the time series. That's why it is not an easy task and many attempts 
are needed before converging towards acceptable solutions. In Fig. [5] we show the result for 
the trading system described above. We observe a correct behavior of the equity, which qual- 
ifies definitely this system. We interpret this as a clear evidence that the time series of the EC 
exhibits features of trend-following, under certain conditions, as encoded in our trading system. 

Finally, in Fig. [61 we check the result on the live-sample period (iii), in 2010-201 1. Here, we 
do not expect any failure, otherwise the full process described above must be rejected. Effec- 
tively, we observe a nice behavior of the cumulative equity (net of fees), much compatible with 
what has been designed on the in-sample. This confirms our statement above on the dynamical 
content of the trading algorithm we have presented here. 

In order to illustrate very simply the gross feature of the model, we present two distributions 
in Fig. |7J We show the trade return spectrum (Fig. [7]- left), in which we recognize a typical 
trend-following system, reminiscent from the standard behavior plotted in Fig. |2] We observe 
also in Fig. [7] (right) that the system is effectively working at high frequency with an average 
duration of trades of 25 minutes. 

From the above discussion, we know that our system is robust against a variation of the 
sampling of the data, for example from 5 minutes to 10 minutes. Also, as the average duration 
of trades is 25 minutes for the nominal system, it makes sense to move the system from 5 
minutes to 10 minutes data sampling and check the results. In order to make the change in an 
optimal way, we have rescaled some parameters such that we have a perfect homotopy between 
the construction at 5 and 10 minutes sampling. Results on the live sample of data (iii) is shown 
in Fig. [8j We observe the good behavior, in accordance (homotopy) with Figj6l as expected. 

A final comment is in order concerning the stress tests and robustness analysis. This study 
ensures statistically the reasonable functioning of the system whatever market regimes and trad- 
ing conditions. We can see the system as an unfolding procedure, transforming the price series 
to a trade series. If we do not control obviously the price series, we have a control on how the 
trade series develop, regardless the behavior of the price series. That's what we have shown 
above. This is a decisive element in the construction of the model. The feed back from the 
cumulative equity itself is an element of the strategy reconstruction, involved with the basic 
conditions of the model as described earlier. Let us add that this is a clear advantage we get 
using this kind of approach for decoding the markets to a certain extend. 
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3 Strategy on the DAX future 



3.1 Data sets and Data treatment 

The idea is to use the same system, as built in section [2l on the DAX future (FDAX) series. 
If we use exactly the same system as defined for Euro future, we obtain a Sharpe ratio of 0.7 
over 10 years of data. This is reasonable but it makes sense to re-optimize some parameters on 
the FDAX series. In the updating process, some conditions are also re-examined and modified 
according to what data (in-sample) requires. We use again a five minutes sampling for the 
FDAX time series, from January 1999 till August 201 1 . This which makes 460k quotes of data. 
We conserve only the close of each quote. As in section [2l this large sample of data points is 
necessary to infer statistical properties with a high confidence level. 

A typical quote of the FDAX is like 5900. The unit of the last digit is what we call a basis 
point. For example, we consider that a price movement from 5900 to 5910 corresponds to a 
price change of 10 basis points. More precisely, if we buy the contract at time Tl (on the quote 
Ql) at 5900 and sell this contract at time T2>T1 (on the quote Q2) at 5910, then this trade 
corresponds to a gain of 10 basis points (without fees). To keep the procedure as close to reality 
as possible we consider fees of two times the slippage, which means that this trade is counted 
in our approach as a trade of 9 basis point (net of fees). 

We follow the analysis process detailed previously. Then, 3 samples of data are defined: 

(i) 1999-2006: in-sample 

(ii) 2007-2009: out-sample 

(iii) 2010-2011: live-sample 

As we start the analysis of the FDAX in 1999, we end up the in-sample in 2006 and not 2007 
as was done for the EC time series. 

3.2 Experimental procedure 

As mentioned above, the basic elements of the algorithm design are the same as the ones used 
for the EC time series in section [2] The main differences concern the treatment of volatility 
ranges. There is a stronger focus on this issue for the FDAX model. 

We have always 8 parameters optimized on the in-sample (i), the system is then validated 
on (ii) and observed on (iii). The optimization is performed in order to achieve the best Sharpe 
ratio on (i). Results are summarized in Fig. [9j We present the behavior with time of the FDAX 
contract itself as well as the cumulative equity of the designed trading system (expressed in 
basis points, net of fees) for the three data samples defined above. 

We observe the nice behavior of the equity, increasing with time, which shows that the 
strategy is profitable and coherent with respect to different market regimes. We observe also 
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that the out- sample (ii) validates the good behavior of the strategy, which is confirmed on the 
recent period 2010-2011. As explained in section [21 we have also guaranteed the robustness 
of the algorithm using a battery of stress tests. In all cases, we must observe that the system 
is stable and robust. This is shown in Fig. [Kjl where we present the Sharpe ratio for all stress 
tests considered. As in section [2l the idea is that the robustness is ensured if we do not observe 
pathological values in the Sharpe ratios, even for the more extreme stress tests. 

We interpret the good result obtained in Fig. [9] as a clear evidence that the time series of the 
FDAX exhibits features of trend-following, under certain other conditions, as encoded in our 
trading system and as already observed on the EC time series. This will be confirmed again on 
the Cacao future in the next section. 

Finally, in Fig. [TH we illustrate explicitly the result on the live-sample period (iii), in 2010- 
2011. Effectively, we identify a nice behavior of the cumulative equity (net of fees), much 
compatible with what has been designed on the in-sample. This confirms our statement above 
on the dynamical content of the trading algorithm we have presented here. 

In order to illustrate very simply the gross feature of the model, we present two distributions 
in Fig. [121 We show the trade return spectrum (Fig. [T2l-left). in which we recognize a typical 
trend-following shape. We observe also in Fig. [121 (right) that the system is effectively working 
at high frequency with an average duration of trades of less that 15 minutes. 



4 Strategy on the Cacao future 

For the Cacao future (CC), we have at our disposal the data series ranging only from 2003 till 
mid-2010. Then, we have chosen to run exactly the EC system on this index, with fees for each 
trade always equal to 2 times the slippage. We do not produce any further optimization. Results 
are presented in Fig. [131 We observe a reasonable equity curve, which proves that the system is 
functioning correctly on this index, uncorrelated to EC and FDAX. 

This confirms the message of this article that the trading system defined above contains 
elements of invariants of financial markets. 



5 Conclusion 

Tools of statistical physics iPTT-PPill have been proven to be efficient in many scientific areas. In 
a similar way for financial time series, knowing that the difference between real markets and 
random walks is very small, a modern statistical multivariate analysis can help to extract this 
difference. This is what is encoded in trading systems. We have shown how to achieve the 
construction of such a system on the Euro future contract at high frequency. A typical element 
of the dynamics of this system is then accessible, namely the trend-following content involved 
in a more complex architecture on volatilities. 

Then, we have produced other examples on completely different markets, largely uncorre- 
lated to the Euro future, the DAX and Cacao future contracts. The same procedure is followed 
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using similar seed ideas and technical inputs. We have shown that similar results can then be 
obtained and we conclude that this is an evidence that some invariants, as encoded in our sys- 
tem, have been identified, on very different markets explored over a 10 years period of time. 
One essential point in our process is that trading models, like the one used in our approach, 
are highly sensitive to non-linear relations in price series. This comes with the multivariate 
data analysis. In this article, we have also focused the discussion on the necessity of a deep 
robustness analysis to ensure the validity of the overall construction. 

An immediate question can be raised concerning the rationale behind this content. Our 
observation is universal in the sense that the same algorithm, for example on EC, is running 
on more that 10 years of data, where the monetary policy has changed several times. Then, 
our approach is not attached to a particular regime of interest rates. There are certainly herding 
behaviors at the origin of the values of parameters encoded in our system. These herding phases 
may appear with strengths governed by certain fear levels, corresponding to volatility domains. 
Also, in some circumstances, nothing special can be said. Finally, a global rationale explanation 
of a given trading system is very complex and probably not unique. This is beyond the scope of 
this article. See ideas in |[T7T - fT9"l . 

In this article, we have completed a pure experimental analysis. The concept of invari- 
ance comes with the observation that we extract similar seed features from largely uncorrected 
financial time series. This is a first step, rooted exclusively on data. 
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Figure 1: Trade return distribution on the Euro future contract using 5 minutes quotes. The 
distributions are illustrated for 2 different years of data, in 2002 and 2005, in order to show the 
similarity in shape. 
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Figure 2: Trade return distribution for a theoretical trend-following system (see text). 
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Figure 3: Behavior with time of the EC contract on the in- sample part (i) of the data series (full 
line) as well as the cumulative equity (dashed line) of the designed trading system (expressed 
in basis points on the right vertical axis, net of fees). The bottom plot corresponds also to 
the running of the trading system, but this time on the randomized in-sample. We observe the 
performance degradation to a well acceptable level. See text for details. 
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Figure 4: Spread of Sharpe ratios corresponding to all stress tests considered in order to ensure 
the robustness of the trading system on the Euro future (see text). 
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Figure 5: Behavior with time of the EC contract on the out- sample part (ii) of the data series (full 
line) as well as the cumulative equity (dashed line) of the designed trading system (expressed 
in basis points on the right vertical axis, net of fees). 




Figure 6: Behavior with time of the EC contract on the live- sample part (iii) of the data se- 
ries (full line) as well as the cumulative equity (dashed line) of the designed trading system 
(expressed in basis points on the right vertical axis, net of fees). 
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Figure 7: Euro future (EC). Left: Trade return distribution for the live-sample (iii). Right: Trade 
duration for the live- sample (iii). 




Figure 8: Behavior with time of the EC contract (sampled in 10 minutes quotes) on the live- 
sample part (iii) of the data series (full line) as well as the cumulative equity (dashed line) of 
the designed trading system (expressed in basis points on the right vertical axis, net of fees). 
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Figure 9: Behavior with time of the FDAX index over the full period 1999-2011 (full line) 
as well as the cumulative equity (dashed line) of the designed trading system (expressed in 
basis points on the right vertical axis, net of fees). We have separated the three data samples 
with dotted lines: (i) 1999-2006 (in-sample), (ii) 2007-2009 (out-sample) and (iii) 2010-2011 
(live-sample). See text for details. 
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Figure 10: Spread of Sharpe ratios corresponding to all stress tests considered in order to ensure 
the robustness of the trading system on the DAX future (see text). 



16 




Figure 11: Behavior with time of the FDAX index on the live-sample part (iii) of the data 
series (full line) as well as the cumulative equity (dashed line) of the designed trading system 
(expressed in basis points on the right vertical axis, net of fees). 
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Figure 12: DAX future. Left: Trade return distribution for the live-sample (iii). Right: Trade 
duration for the live- sample (iii). 



17 




Figure 13: Behavior with time of the CC index over the period 2033-2010 (full line) as well as 
the cumulative equity (dashed line) of the designed trading system (expressed in basis points on 
the right vertical axis, net of fees). 
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